Explosion of a Collapsing Bose-Einstein Condensate 
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We show that elastic collisions between atoms in a Bose-Einstein condensate with attractive 
interactions lead to an explosion that ejects a large fraction of the collapsing condensate. We study 
variationally the dynamics of this explosion and find excellent agreement with recent experiments 
on magnetically trapped 85 Rb. We also determine the energy and angular distribution of the ejected 
atoms during the collapse. 
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Introduction. — Ever since the achievement of Bose- 
Einstein condensation in an atomic gas with attractive 
interactions JjJ, it has been an important objective to 
study in detail the condensate collapse that is predicted 
to occur HH if the number of condensate atoms exceeds 
a certain maximum number determined by the strength 
of the attractive interactions In the pioneering ex- 
periments of Bradley et al. this objective was very diffi- 
cult to achieve for two reasons. First, the use of doubly 
spin-polarized atomic 7 Li results for these experiments in 
a maximum number of condensate atoms of only about 
1400 atoms, which is too small to allow for nondestruc- 
tive imaging of the collapse dynamics. Second, the in- 
herent stochastic nature of the growth and collapse cy- 
cles that occur because one aims at evaporatively cooling 
far below the critical temperature || , prevents sufficient 
control over the initial conditions of the condensate to 
perform a sequence of destructive measurements. How- 
ever, a statistical analysis of the condensate growth and 
collapse cycles has nevertheless revealed important infor- 
mation on the collapse process || , and new experimental 
techniques are presently being applied to overcome in 
particular the latter of these problems . 

Complementary to these developments, Cornish et al. 
have recently been able to overcome both problems in 
an ingenious experiment with spin-polarized atomic 85 Rb 
H . In this experiment one makes use of the fact that the 
/ = 2, 771/ = —2 hyperfme ground state of 85 Rb has a 
so-called Feshbach resonance |J at a magnetic field of 
about 1.55 mT, which offers the opportunity to mag- 
netically tune the interatomic interactions from strongly 
repulsive to very attractive. As a consequence Cornish 
et al. were able to first produce a large, stable and es- 
sentially pure condensate, and then suddenly switch the 
interactions from repulsive to attractive to induce a col- 
lapse and observe its properties. 

Perhaps the most surprising outcome of this experi- 
ment is that during the collapse an explosion occurs that 
ejects a large number of relatively highly energetic atoms 
from the condensate. A possible mechanism that imme- 
diately comes to mind for the production of these en- 
ergetic atoms are inelastic collisions between condensate 



atoms that flip a spin and, therefore, convert Zeeman 
energy into kinetic energy. However, a simple estimate 
shows that this mechanism generally leads to particles 
with much too high energies and is also incapable of 
ejecting so many atoms from the condensate. Similarly 
it appears that elastic collisions between a condensate 
atom and an atom in the thermal cloud, which in prin- 
ciple cannot experimentally be excluded to be present, 
also occur too infrequent to explain the experimental re- 
sults. We are thus faced with the important theoretical 
task to determine what physical mechanism is responsi- 
ble for the observed explosion. It is the main purpose 
of this Letter to point out that elastic collisions between 
two condensate atoms can provide such a mechanism. 

Elastic Condensate Collisions. — In detail, the scat- 
tering process that we have in mind is a collision between 
two condensate atoms by which one of the atoms is stim- 
ulated back into the condensate again but the other atom 
is ejected out of the condensate. It is important to realize 
that such a process is forbidden for a homogeneous con- 
densate due to momentum conservation. However, for a 
trapped gas this is no longer true, because the condensate 
now occupies a band of low- lying momentum states. Note 
also that momentum and energy conservation ensure that 
the ejected atoms will automatically have, as compared 
to inelastic spin- flip collisions, small kinetic energies and 
can thus remain trapped as seen in experiment. To see if 
the above mechanism can also explain the large amount 
of ejected atoms, we need to calculate the rate associated 
with this process. 

This can be most easily achieved by using the fact that 
we are dealing with a two-body interaction process. As a 
result the rate of change of the total number of conden- 
sate atoms N c (t) is related to the instantaneous interac- 
tion energy of the condensate E 1 ^ (i) by 



dN c (t) 
dt 



flm^t)] 



(1) 



In the Bogoliubov theory of the dilute Bose gas we simply 
have that 
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which in momentum space reads 

x(27r) 3 ,5(k4+k3-k 2 -k 1 )r 2B 0(k a ,t)^(ki,t) . (3) 

Here the two-body T(ransition) matrix equals T 2B = 
An aft 2 /m in terms of the s-wave scattering length a and 
the mass m of the 85 Rb atoms, and 4>(x, t) is the conden- 
sate wave function. In this approximation the condensate 
interaction energy is thus purely real and the number of 
condensate atoms is conserved. 



However, a microscopic derivation jlO|,[ll| shows that 
in Eq. (g) the expression (27r) 3 <5(k 4 + k 3 - k 2 - ki)T 2B 
should in principle be replaced by the many-body T- 
matrix element T MB (k4, k3, k 2 , ki), whose real part is at 
low temperatures indeed well approximated by the Bo- 
goliubov result (27r) 3 <5(k 4 + k 3 - k 2 - ki)T 2B . However, 
the many-body T matrix describes also real (incoherent) 
collisions taking place in the gas and as a result acquires 
an imaginary part that can be obtained from an optical 
theorem [|l0| . Being interested in the ejection of relatively 
highly energetic particles, we can in first instance neglect 
the mean-field effects on the corresponding intermediate 
state to obtain 



Im[T MB (k 4 ,k 3 ,k 2 ,k 1 )] 

= -2n(T 2B ) 2 / -—3 <5(e(k 2 )+e(k 1 )-e(k)-e(k 2 +k 1 -k)) 0(k* + k 3 - k, t)0*(k a + k x - k, t) , (4) 



where e(kj) = /i 2 k 2 /2m, Kk is the momentum of the 
ejected atom, and the dependence on the condensate 
wave function is caused by the effect of Bose enhance- 
ment of scattering into already occupied states. Interest- 
ingly, the latter is here nondiagonal in momentum space 
due to the inhomogeneity of the condensate. 

Notice that from a field-theoretical point of view, we 
have in the above manner just arrived at an evaluation of 
the Feynman diagram drawn in Fig. 0. We thus conclude 
that we are effectively dealing with a three-body process. 
Indeed, by taking the functional derivative 8/8<fi*(x., t) of 
our result for the condensate interaction energy, we can 
obtain an imaginary three-body correction term to the 
Gross-Pitaevskii equation of the condensate, which pro- 
vides another way of deriving the desired ejection rate. 
It is also worth mentioning that we qualitatively expect 
mean-field effects to reduce the above ejection rate in the 
case of a gas with positive scattering length, because it 
then costs energy to remove a particle from the conden- 
sate. Such an effect, however, does not occur for the 
negative scattering length case of interest here. 



Collapse Dynamics. — To make further progress we 
need to determine the condensate wave function during 
the collapse. We do this by solving the Gross-Pitaevskii 
equation variationally, taking for the wave function an 
anisotropic Gaussian with time-dependent widths de- 
noted by q x (t), Qy(t), and q z {t). In detail we use, apart 
from a for our purposes irrelevant phase factor, 



(x,i) = v/iV, 



x exp ■ 



n 



1 



1/4 



T9?(t) 



. mqj (t) dqj (t) ' 



dt 



(5) 



The reason for this Gaussian ansatz is that it is the appro- 
priate description of the condensate after the explosion 
has occurred and the remnant contains only a relatively 
small amount of atoms P,p"2"|. Moreover, it is well-known 
that a Gaussian ansatz gives excellent results for the col- 
lective modes of the condensate, even when the number 
of atoms in the condensate is so large that the Gaus- 
sian approximation for the condensate wave function is 
strictly speaking no longer valid p3| . We thus expect it 
to give physically sensible results for the dynamics of the 
collapse at all times. 

Within this approach the explosion of the collapsing 
condensate is described by the 'classical' equations of mo- 
tion 



V 



FIG. 1. The imaginary part of this Feynman diagram de- 
termines the rate for ejection of atoms out of the condensate 
due to elastic condensate collisions. The dashed lines corre- 
spond to condensate atoms and the full line to a nonconden- 
sate atom. 



with a potential energy that is equal to 



V( q ;N c )=J2 



2 2 



2mq 2 



2 ah 2 N c 
7T mq x q y q z 
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for an anisotropic harmonic trapping potential with 
spring constants mux. Moreover, these equations are 
coupled to the rate equation for the number of conden- 
sate atoms that follows from Eqs. ([!]) and (Q). In the 
anisotropic case this cannot be worked out fully analyti- 
cally and we have to determine the ejection rate numeri- 
cally. However, in the isotropic case we simply find that 



dN c {t) _ 2V5a 2 hN^{t) 
dt ~ ~ mq 4 (t) 



(8) 



The dependence of the right-hand side on N c (t) and 
q(t) is somewhat unusual, because if we had discussed 
the effect of inelastic two-body processes the rate of 
change of the number of condensate atoms would be 
proportional to N% (t) / q 3 (t) . The additional factor of 
N c (t) is easily understood and reflects the fact that 
the ejection of the atoms is Bose stimulated. As men- 
tioned previously, our mechanism therefore effectively 
behaves as a three-body process. The additional factor 
of l/q(t) is more subtle and shows that, if the con- 
densate collapses and the wave function becomes more 
spread out in momentum space, more atoms can satisfy 
the energy and momentum constraints for the ejection. 
Having said this, we should make sure that the domi- 
nant contribution to the integration in the right-hand 
side of Eq. (|J) comes from momenta h\i that lie out- 
side the band of momenta occupied by the condensate. 
We can most conveniently achieve this by including in 
the integrand the factor (1 - |0(k, £)| 2 /|0(O, i)| 2 ), which 
smoothly interpolates between for k • q <C 1 and 1 
for k ■ q > 1. In the following we always use this 
smooth cut-off, because atoms with momenta below the 
cut-off are most likely stimulated back into the con- 
densate and do not contribute to the ejection rate. It 
has, however, only a small quantitative effect as can be 
seen explicitly in the isotropic case, where the right-hand 
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FIG. 2. The radial and axial widths of the condensate dur- 
ing a collapse induced by a sudden change of sign in the scat- 
tering length. At the origin of the time axis the scattering 
length vanishes. See text for more details. 
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FIG. 3. The number of atoms in the condensate during a 
collapse induced by a sudden change of sign in the scattering 
length. The inset shows the ejection rate. 



side of Eq. ( p[) j ust becomes multiplied by the factor 
(1 - (8/25)^/2/5) ~ 0.8. 

We have solved the above set of coupled equations for 
the conditions of the experiment performed by Cornish et 
al. |@| . In this experiment one uses a cigar-shaped mag- 
netic trap with a radial frequency of uj r = lu x = uj y — 
2ir x 17.5 Hz and an axial frequency of oj z — 2tt x 6.8 
Hz. Moreover, one first makes a large condensate with 
about 4000 atoms and a positive scattering length of 2500 
do. The scattering length is then within 0.5 ms changed 
to —60 ao by means of a linear ramp in the magnetic 
bias field. The outcome of our simulation of this exper- 
iment is summarized in Figs. |^ and |^, and appears to 
be in excellent agreement with the preliminary experi- 
mental data |14|. From these figures we see that the 



condensate collapses first in the radial direction in ap- 
proximately 7r/(2cj r ) ~ 14 ms. During the last part of 
this collapse an explosion occurs in which about 3/4 of 
the initial number of atoms is expelled from the conden- 
sate. As a result the number of condensate atoms is now 
less than the maximum number of atoms possible to have 
a metastable condensate, and the radial collapse of the 
condensate is turned into a large amplitude oscillation. 
In principle we can discern at each inner turning point of 
the radial and axial oscillation an increased loss of con- 
densate atoms, but these are never as dramatic as in the 
first time, since the condensate is now no longer unstable. 

Explosion Process. — To explore the physics of the ex- 
plosion further, we have also determined the energy and 
angular distribution of the ejected atoms, which have not 
yet been examined in detail experimentally. In Fig. |] we 
show the energy distribution at two times during the first 
radial collapse of the condensate. What is most striking 
is that at later times the distribution is much broader. 
This is again caused by the condensate wave function 
being much more spread out in momentum space at the 
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FIG. 4. The energy distribution of the atoms ejected from 
the condensate at two times during the first radial collapse. 
We have used a dimensionless energy variable by scaling the 
energy with hcu = 7i(tj^aj z ) 1 ' /3 ~ 0.61 nK. 

later time, so that the ejected atoms can gain much more 
energy from the condensate collisions. In Fig. || we have 
also depicted at various times the angular probability dis- 
tribution, which turns out to be essentially independent 
of the energy of the ejected atoms but is clearly seen 
to depend strongly on the ratio of the radial and axial 
widths of the condensate. In particular, at the peak of 
the explosion the distribution is very anisotropic and al- 
most no atoms are being ejected along the z-axis. This 
can be understood from the fact that the condensate has 
at that time a very elongated cigar shape. As a result 
both the total momentum of the two colliding conden- 
sate atoms as well as the momentum of the atom that is 
stimulated back into the condensate, always have to be 
directed almost perpendicular to the z-axis. The same is 
therefore true for the momentum of the ejected atom. 
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FIG. 5. Polar plot of the angular probability distribution 
of atoms ejected from the condensate at three different times 
but at the same energy. 

Discussion. — Although we have focused in this Letter 
on the recent exciting experiments with 85 Rb, we believe 



that the new condensate loss mechanism that we have 
proposed here is also important for the experiments with 
7 Li and may resolve the existing discrepancy between 
theory and experiment in that case |(| . In this respect it 
should be noted that the collapse occurring in the 85 Rb 
experiments is physically quite different and in a sense 
not so violent than the one occurring in 7 Li. Roughly 
speaking the difference is that in the former case the large 
condensate ejects so many particles that the remnant of 
the explosion corresponds to a metastable condensate, 
whereas in the latter case it is precisely this metastable 
condensate that collapses to even smaller sizes. Possibly 
another important difference is that in the experiments of 
Cornish et al. no thermal component is visible, whereas 
the experiments of Sackett et al. are close to the critical 
temperature and, therefore, a large thermal cloud is con- 
stantly feeding the condensate. We intend to come back 
to a detailed theoretical treatment of these interesting 
issues in a future publication. 
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